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A method providing real-time determination of a thermal profile of the drilling 
fluid in a well from three measuring points available in the field, i.e. the injection 
temperature, the outlet temperature and the bottomhole temperature is 
disclosed. The form of the profile between these three points is defined by a type 
curve representative of the thermal profiles in a well under drilling, estimated 
from physical considerations on thermal transfers in the well. 



ABSTRACT 




BACKGROUND OF THE INVENTION 

Field of the Invention 

[0001] The present invention relates to a method of determining the 
thermal profile of a drilling fluid in a well. 

[0002] During drilling, which is the mud injected into the drill string of the 
well and which flows back through the corresponding annulus, undergoes 
great temperature variations. The fluid can encounter temperatures that can 
range from 2^*0 for deep offshore wells to more than 180°C for very hot 
wells. Many mud properties, such as rheology or density, depend on the 
temperature. Calculation of the pressure losses during drilling can therefore 
be improved if an estimation of the temperature profile in the well is known. 
It is therefore important to be able to predict the temperature profile in the 
flowing mud from well data and mud characteristics. 

[0003] Measurement of the thermal profile of the fluid in a well under 
drilling would require complete instrumentation of the well, that is installation 
of evenly spaced out detectors in the drill string and in the annulus. allowing 
temperature measurement at various depths. However, installing such a 
measuring system entails too many constraints ; only localized 
measurements picked up by devices mounted in the drill string allow 
knowing certain temperature points on the path of the drilling fluid. 



[0004] In the face of this lack of data, analytic models based on heat 
transfer equations have been developed to evaluate the thermal profiles of 
the fluid along the well under drilling. Some of these analytic models are 
implemented in softwares and allow providing an estimation of thermal 
profiles from a data which may be difficult to obtain. Thus, knowing the 
characteristics of the site and of the drilling equipment, by giving a value of 
the temperature of the fluid at the well inlet, these softwares can predict the 
temperature profile of the drilling fluid. 

[0005] However, a comparison between the results obtained with analytic 
methods and the measurements obtained in the field shows that there can 
be great differences. Furthermore, the complexity of softwares using 
numerical calculation methods makes real-fime implementafion thereof 
difficult. 

[0006] On the other hand, a study of the bibliography on thermal models 
shows a similarity in the form of temperature profiles in most cases, which 
turns on three points : inlet temperature, outlet temperature and bottomhole 
temperature. 

[0007] The aim of this study is thus to provide a method allowing real- 
time determination of a thermal profile in the mud from three measuring 
points available in the field, that is the injection temperature, the outlet 
temperature and the bottomhole temperature measured by a detector 
mounted on the string. The form of the profile between these three points is 
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represented by a type curve representative of the thernnal profiles in a well 
under drilling, estimated from physical considerations on thermal transfers in 
the well. 

SUMMARY OF THE INVENTION 

[0008] The method of determining the thermal profile of a drilling fluid 
circulating in a well under drilling according to the invention is defined by the 
successive stages as follows : 

a) determining a general expression 61 for the thermal profile of the fluid 
within the drill string in the well and a general expression 92 for a thermal 
profile of the fluid in the corresponding annulus. using the heat propagation 
equation that takes into account a thermal profile of the medium surrounding 
the well, 

b) measuring the temperature of the fluid at the well inlet. T1, at the well 
bottom, T2, and at the well outlet, T3, 

c) modifying the expressions 91 and 92 to meet the temperature boundary 
conditions T1 , T2 and T3. 

d) detennining the thermal profile of the drilling fluid as a function of the 
depth. 

[0009] In order to obtain, in real time, a temperature profile with the 
method presented above, stages b), c) and d) can be repeated. 

[0010] According to the method of the invention, in stage a), general 
expressions 91 and 92 can comprise unknown constants, and in stage c), it 
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can be determined that expressions 91 and 92 meet the temperature 
boundary conditions T1 , T2 and T3 by determining the unl<nown constants. 

[0011] In order to determine a general expression 91 for the thermal 
profile of the fluid within the drill string in the well and a general expression 
92 for a thermal profile of the fluid in the corresponding annulus, it is 
possible, according to the method of the invention, to use in stage a) the 
heat propagation equation that takes into account at least the thermal 
equation of the medium surrounding the well, the flow rate of the fluid and 
the balance of the thermal exchanges undergone by the fluid, the thermal 
exchanges comprising at least the exchanges between the ascending and 
descending drilling fluid, and/or to use the equation of heat propagation in a 
homogeneous medium on a cylinder of infinite height centered on the well, 
the cylinder comprising the drill string that guides the descending fluid and 
the annulus around said drill string, which guides the ascending fluid. 

[0012] According to the method of the invention, general expressions 91 
and 92 obtained in stage a) can be split up into several independent 
equations and, in stage c), it can furthermore be determined that the profiles 
and the derivatives of the thermal profiles of the fluid within the drill string 
and in the corresponding annulus are continuous. 

[0013] The method according to the invention can notably be used to 
calculate the pressure drops of the drilling fluid circulating in a well under 
drilling, or in another application, to determine the zones of hydrate 
formation in the fluid during drilling. 
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[0014] In relation to the methods for determining the thermal profile of a 
drilling fluid in a well according to the prior art, the present invention notably 
affords the following advantages : 

- the temperature profile obtained is more accurate because it is 
determined from three drilling fluid temperature measurement points while 
keeping an analytic expression of the thermal profile between the measuring 
points which is physically justified, 

- since the temperature measurements are performed all the time, the 
method allows obtaining the temperature profile in real time and to observe 
the evolution thereof with time. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] Other features and advantages of the present invention will be 
clear from reading the description hereafter of non-limitative examples, with 
reference to the accompanying drawings wherein : 

- Fig. 1 diagrammatically shows the architecture of a well under drilling, 

- Figs. 2, 3 and 4 show the form of the temperature profile of the drilling 
fluid in a vertical onshore well, 

- Fig. 5 shows the form of the temperature profile of the drilling fluid in a 
vertical offshore well, 

- Fig. 6 shows the form of the temperature profile of the drilling fluid in a 
deviated offshore well. 
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- Fig. 7 shows the evolution as a function of time of the temperature profile 
of the drilling fluid in a vertical offshore well. 

DETAILED DESCRIPTRION OF THE PREFERRED EBMODIMENTS 

[0016] It is possible to give an analytic expression for the thermal profile 
in the well and the annulus by using quite simple heat exchange 
considerations, that is the heat propagation equation. 

[0017] This model is based on the establishment of the heat balances in 
the well. According to a first approach, only the steady states are 
considered (the drilling mud flow is assumed to be stabilized for some time 
so that the temperatures no longer evolve). Certain hypotheses are 
necessary for calculation : the heat exchanges are measured in a plane 
perpendicular to the laminar flow of the mud, the various constants are 
assumed to be independent of the temperature, and finally the influence of 
the temperature of the medium surrounding the well shows on an apriori 
selected useful diameter Rf. 

[0018] It is then sufficient to use the heat propagation equation in a 
homogeneous medium on a cylinder of infinite height centered on the well 
shown in Figure 1. In each well section, the heat loss equality is written by 
considering two temperature functions : ei(z) within the drill string and Qziz) 
in the annulus. 
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[0019] Let 

Gf be the temperature of the formation, 

Xf the thermal conductivity of the medium surrounding the well, 

Xa the thermal conductivity of the tubing (metal), 

Cp the heat-capacity rate of the drilling fluid, 

R1 the inside radius of the drill string. 

R2 the outside radius of the drill string, 

Rt the radius of the annulus, 

Rf the effective radius (for heat supply) around the well, 
D the flow rate of the drilling fluid, 
p the density of the drilling fluid. 

[0020] The heat balances per unit of depth are as follows : 

- Heat supplied by the medium surrounding the well to the fluid in the 



in A J. 






In 







annulus : 

- Heat carried from the fluid in the annulus to the fluid within the drill 
string : 

- Heat accumulated by the fluid in the drill string and in the annulus : 



R 



1 J 



7 



Qt = -D.p.CpAGi 
Qa = D.p.CpAe2. 
[0021] The heat balances lead to the following system 

Qt = Q2 

Qa = Qi + Q2. 

That Is 



ie. 



dz 
dz 
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[0022] These equations are solved by diagonallzatlon and matrix 
inversion , and they lead to the following results : 



a 



with 



A = 



DpC, In 
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Gf = a.z+eo being the thermal equation of the medium surrounding the well 
and a the thermal gradient. 

Ki and K2 are the integration constants depending on the boundary 
conditions. 

[0023] It is thus possible, by means of some simplifying hypotheses, to 
obtain an analytic expression for the temperature profile of the drilling fluid 
in a well. If all the parameters are known, by giving the inlet temperature 
and by writing that temperatures 01 and 82 are equal at the well bottom, the 
profile is entirely determined. The most commonly used softwares use this 
type of predictive procedure. However, a study of the results of the models 
compared to field data shows how difficult it is to use these estimations in a 
predictive way. 

[0024] In the present invention, the system is based on the knowledge of 
three measuring points in the field : inlet temperature, outlet temperature 
and bottomhole temperature. In order to estimate the thermal profile in the 
well from the three measurements consisting of the surface injection and 
outlet temperatures and the bottomhole temperature (inside or outside the 
drill string), the method according to the invention connects the three 
measuring points by a general expression representative of the evolution of 
a thermal profile in a wellbore, as obtained according to the method 
described above. 

[0025] Therefore the equations obtained by means of heat exchange 
calculations as follows are used : 
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9, (z) = -K^Be'" - K^Be'^' +9^-- 

B 

9,iz) = -K,{B + r,)e'" - K^{B + r^)e'^' + 9^ 

[0026] According to the invention, these curve fornns are adjusted to the 

0,(z) - -K,Be''' -K^Be'^' +9^- — 

B 

0^(z) = -K,(B + r,)e'' ' - K,{B + r^K'" + 

three measuring points of the driHing fluid temperature at the inlet, T1 , at the 
well bottom, T2, and at the well outlet, T3. In order to use these three 
measuring points as boundary conditions, the two equations (in the drill 
string and in the annulus) are decoupled by using different integration 
constants while keeping the general expression. Two general expressions 
of the temperature profile are obtained in the drill string, 91, and in the 
annulus, 92, which have a physical significance but which comprise two 
degrees of freedom. Thus, expressions 91 and 92 can be adjusted by filing 
the degrees of freedom in order to meet the temperature conditions T1 , T2 
and T3. Therefore the equations in the pipes and in the annulus have the 
form as follows : 

[0027] Four integration constants Ki, K2. K3 and K4 are obtained rather 
than two, which requires four boundary conditions to determine the 
temperature profile. These four boundary conditions are : measurements of 
the well inlet, bottomhole and outlet temperature, and an equality condition 
at the bottom between the temperature in the drill string. 91. and the 
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temperature in the annulus. 92. The profile can be adjusted at any time to 
pass through the measuring points : therefore a real-time estimation of the 
thermal profile is obtained. Programming by means of a spreadsheet type 
software allows ready obtaining of the real-time evolution of the profile. 

[0028] Figs. 2, 3 and 4 respectively show the temperature profile of the 
drilling fluid in a vertical onshore well at a flow rate of 500 l/min, 1000 l/min 
and 2000 l/min. The analytic expression determined allows simple 
calculation of the temperature T in Celsius degrees of the fluid in the drill 
string (curve 61) and in the annulus (curve 02) as a function of the depth P 
in meter. The analytic expression depends on several parameters that can 
be fixed from the start. By default typical values of these parameters are 
used. To determine the temperature profile of Figs. 2, 3 and 4, the 
geothermal gradient a is assumed to be constant to correspond to the 
onshore situation of the well. The temperature profile is entirely determined 
by measuring the temperature, 20°C at the inlet. 35°C at the bottom and 
24°C at the outlet of the well. 

[0029] The case of the vertical offshore well can be addressed by 
considering that the geothermal profile of the medium surrounding the well 
is divided in two domains : let Om be the thermal profile of the sea and Gs 
the thermal profile of the ground. Thermal gradient a is assumed to be 
constant in each domain, but discontinuous exist from one domain to the 
other. Let am be the thermal gradient of the sea and as the thermal gradient 
of the ground. Two series of equations (one for each domain) are 
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considered for each general expression in the pipes and in the annulus. 
Four decoupled equations are obtained which represent the thermal profile 
of the drilling fluid in the well. Equation 911 (z) corresponds to the 
temperature profile in the drill string in the sea, ei2(z) corresponds to the 
temperature profile in the drill string in the ground, 921 (z) corresponds to the 
temperature profile in the annulus in the ground and 922(z) corresponds to 
the temperature profile in the annulus in the sea, 91 1 being independent of 
912 and 921 being independent of 922 : 

CCS 

^„(z) = -K,Be'' ' -K,Be'''+0„- — 

0„(z) = -K,{B + r,)e'' ' - K,iB + r,)e''' + 0, 
0,,{z) = -K,{B + r,)e^'-^ -K,{B + r,)e''-' + 0„ 

[0030] This brings the number of integration constants to eight (Ki to Ks). 
The boundary conditions are then : measurements of the inlet, outlet and 
bottomhole temperature, equality condition at the bottom between the pipe 
temperature and the annulus temperature, to which the continuity of the 
thermal profiles is added in the drill string and in the annulus at the junction 
of the two domains and the continuity of the derivative of the thermal profiles 
in the drill string and in the annulus at the junction of the two domains. 
Similarly, it is then possible to obtain in real time a physically realistic 
thermal profile passing through the measuring points. Fig. 5 shows the 
thermal temperature profile of a drilling fluid in an offshore well from the four 
equations 911, 912, 921 and 922. The fluid circulates at 500 l/min and the 
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temperatures measured are 20-0 at the inlet. 15°C at the outlet and 30"C at 
the well bottom. The thermal gradients are selected constant In each 
domain crossed by the well. 

[0031] Deviated wells represent the majority of the current wellbores. 
The physical problem is not fundamentally different and it can be handled in 
the same way as offshore wellbores : the well just has to be divided into two 
domains, each domain being characterized by a different thermal gradient 
corresponding to the medium surrounding the well. In the case of a deviated 
well, the depth corresponds to the distance covered along the well 
trajectory. General expressions 91 and 62 representative of the thermal 
profile are each split up into two independent equations. The vertical part is 
characterized by the thermal gradient a of the medium surrounding the well, 
the deviated part is characterized by an equation of the thermal profile of the 
medium surrounding the well Qa = a.sin((|)).z+eo, <t) being the angle of 
inclination. The same boundary conditions (measurements of the 
temperature at the inlet, at the outlet and at the bottom of the well, equality 
at the bottom between the pipe temperature and the annulus temperature, 
and continuity of the thermal profiles and of the derivative of the thermal 
profiles in the drill string and in the annulus at the junction of the two 
domains) then allow solving the equations and obtaining the expression of 
the temperature profile in the pipes and in the annulus. 

[0032] It is possible to combine the procedure applied for a vertical 
offshore well and the procedure applied for a deviated onshore well in order 
to determine the temperature profile in an offshore well when the direction of 
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the hole is deviated in the ground. The domain is divided into three different 
domains : let Gm be the thermal profile of the vertical domain in the sea, Gs 
the thermal profile of the vertical domain in the ground and Gd the thermal 
profile of the deviated domain in the ground. Fig. 6 shows the thermal profile 
in a deviated offshore well. The fluid circulates at 500 l/min and the 
temperatures measured are 20°C at the inlet, 23°C at the bottom and 15°C 
at the outlet of the well. 

[0033] According to the same method as that used for the vertical 
offshore well or for the deviated onshore well, it is possible to determine the 
thermal profile of a vertical onshore well whose formation thermal gradient 
changes as a function of the depth. The well is divided into domains 
characterized by a thermal equation of the medium surrounding the well. 
General expressions G1 and G2 representative of the thermal profile are then 
each split up into as many independent equations as there are different 
domains. The same boundary conditions (measurements of the temperature 
at the inlet, at the outlet and at the bottom of the well, equality at the bottom 
between the pipe temperature and the annulus temperature, and continuity 
of the thermal profiles and of the derivative of the thennal profiles in the drill 
string and in the annulus at the junction of the two domains) then allow 
solving the equations and obtaining the expression of the temperature 
profile in the pipes and in the annulus. 

[0034] By repeating the calculation allowing obtaining the expression of 
the temperature profile of the drilling fluid upon each new temperature 
measurement, a representation of the temperature profile evolving with time 
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is obtained. Fig. 7 shows the evolution of the temperature profile of the 
drilling fluid in an offshore well in the course of time. The graph in the upper 
part of Fig. 7 shows the evolution as a function of time t in seconds of the 
flow rate parameter D in l/min of the drilling fluid, and of the temperature 
parameter T in °C of the drilling fluid at the inlet, T1 , at the bottom, 12, and 
at the outlet, T3, of the well. The three graphs in the lower part show the 
temperature profile at three different times and allow to observe the 
evolution thereof. 

[00351 Knowledge of the thermal profile of the drilling fluid at any time 
allows real-time calculation of the pressure drops in the well by taking into 
account the thermal effects. This gives a better estimation of the bottomhole 
pressures and of the injection pressures for complex wells. 

[0036] Another use of real-time determination of the thermal profile of the 
drilling fluid is hydrate formation prevention. Hydrates fonm under low 
temperature and high pressure conditions, conditions which are met notably 
in deep offshore wells at the ground/sea interface. Knowledge of the 
temperature profile allows determination of the zones where the 
temperature of the drilling fluid is below the minimum value from which 
hydrates form, then to react accordingly, for example by raising the flow rate 
or by heating the fluid in order to prevent this formation of hydrates. 
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